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In the framework of constituent quark model, the effect of hidden color channels on the nucleon- 
nucleon (AW) interaction is studied. By adjusting the color confinement strength between the 
hidden color channels and color singlet channels and/or between the hidden color channels and 
hidden color channels, the experimental data of S to I partial-wave phase shifts of AW scattering 
can be fitted well. The results show that the hidden color channel coupling might be important 
in producing the intermediate-range attraction of AW interaction. The deuteron properties and 
dibaryon candidates have also been studied with this model . 
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I. INTRODUCTION 

The study of nucleon-nucleon (A/A/) interaction has 
lasted over seventy years. The quantitative description 
of AW interaction has been achieved in the one-boson- 
exchange (OBE) models, the chiral perturbation the- 
ory(ChPT) and quark models. The x 2 /dof~ 1 for more 
than 2000 data has been obtained in meson exchange 
model [H, 0] and <2 in quark model 0. 

In the OBE model the long-range part of the 
AW interaction is attributed to one-pion-exchange. The 
short-range part is described by p, w-meson exchange or 
phenomenological repulsive core. While the a-meson 
exchange is responsible for the intermediate-range at- 
traction. Phenomenological form factors are needed to 
achieve the quantitative description of the NN interac- 
tion data. In the chiral perturbation theory Q , the multi- 
7r's are exchanged between two nucleons. The short range 
part related to the nucleon internal structure is modeled 
by the contact terms with phenomenological low energy 
constants. The theory can give a quantitative description 
of the low-energy NN scattering below the 7r production 
threshold. It is hard to extend this model to higher en- 
ergy, the very interesting resonance region of NN scat- 
tering. 

With the advent of quantum chromodynamics (QCD), 
it is expected to describe the NN interaction from the 
fundamental degree of freedom of QCD, quark and gluon. 
Recently, lattice QCD calculation has achieved a quali- 
tative description of NN interaction However it is 
still far from the quantitative description. The QCD- 
inspired quark models are useful in describing the NN 
interaction with the fundamental quark-gluon degree of 
freedom. The most popular and successful one is the con- 
stituent quark model. Where the non-perturbative (color 
confinement and spontaneous chiral symmetry breaking) 
and perturbative properties of QCD are incorporated 
into the model by introducing the phenomenological con- 
finement potential, Goldstone-boson exchange and effec- 
tive one gluon exchange between the massive constituent 
quarks 5] . Almost in all realistic quark models aimed to 



describe the AW interaction, the short-range repulsion 
of AW interaction is described by one-gluon-exchange 
and quark anti-symmetrization. The long-range part is 
described by 7r-meson exchange which is the same as 
the OBE and chiral perturbation theory approaches. To 
describe the intermediate-range part, the a-meson ex- 
change is employed in most quark model approaches. 
The only one exception is the quark delocalization color 
screening model (QDCSM). Where the quark delocaliza- 
tion and color screening effect between interacting quarks 
within different quark clusters are employed [6| to de- 
scribe the intermediate range attraction which is sim- 
ilar to the molecular covalent bond. To develop such 
a molecular covalent bond like model is because of the 
outstanding fact that the molecular force and nuclear 
force are similar except the energy and length scale dif- 
ference 0] ■ Also because of the existence of a meson is 
not sure for long. Recently BES collaboration reported 
the observation of cr-meson, which is appeared as 7T7t S- 
wave resonance Q. However, the calculation of the corre- 
lated 7T7T exchange between two nucleons can not obtain 
enough attraction Q as the phenomenological tr meson 
exchange did. The recent QDCSM calculation, on the 
other hand, showed that the quark delocalization and 
color screening mechanism is quantitatively equivalent 
to the phenomenological a meson exchange in describ- 
ing the AW intermediate range attraction jlO]. In ChPT 
there is also no a meson exchange. In addition, by intro- 
ducing the multi-body color confinement interaction [TTI j , 
or by incorporating the hidden color channels in the cal- 
culation [12j, the intermediate-range attraction can also 
be obtained to some extent. Therefore the mechanism 
of the AW intermediate-range attraction is still an open 
question. 

In this work, an alternative approach for AW inter- 
action is studied. The hidden color channels ignored in 
the prevailing quark model calculations of AW interac- 
tion is included. Accordingly the confinement potential 
between different channels is modified as follows: the or- 
dinary confinement is used for the quark-pairs within the 
same nucleon and the color singlet channels whereas a 
multiplying factor is introduced for the confinement po- 
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tential between the quark pairs if a hidden color channel 
is involved. The aim is to test if the color screening phe- 
nomenology used in QDCSM is an effective description 
of the hidden color channel coupling. The details of this 
model approach will be explained in next section. The 
NN scattering phase shifts obtained in this approach are 
confronted with experimental data and compared with 
ChQM and QDCSM approaches. The equivalence of 
these three quark models in describing the NN scattering 
data has been confirmed. The deuteron properties and 
dibaryon candidates are also studied with this model. 

The structure of this paper is as follows. A brief intro- 
duction of three quark models used is given in section II. 
Section III devotes to the numerical results and discus- 



sions. The summary is shown in the last section. 



II. THREE QUARK MODELS 

A. Chiral quark model 

The Salamanca version of ChQM is chosen as the rep- 
resentative of the chiral quark models. It has been suc- 
cessfully applied to hadron spectroscopy and NN inter- 
action. The model details can be found in Ref. 13]. Only 
the Hamiltonian and parameters are given here. The 
ChQM Hamiltonian in the nucleon-nucleon sector is 
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Where is quark tensor operator, Y(x), H{x) and G(x) 
are standard Yukawa functions T c is the kinetic en- 
ergy of the center of mass, a c h is the chiral coupling con- 
stant, determined as usual from the 7r-nucleon coupling 
constant. All other symbols have their usual meanings. 
The parameters of this ChQM Hamiltonian are given in 
Table H 



B. Quark delocalization color screening model 

The model and its extension were discussed in detail in 

Ref. mm. 

Its Hamiltonian has the same form as Eq.(l), 



but without a meson exchange and a phenomenological 
color screening confinement potential is used, 

rCXS 



-a c X i -X j [f(r ij ) + Vo] + Vl, 



if i, j occur in the same 
baryon orbit, 

if i, j occur in different 
baryon orbits. 



(2) 



Here, fi is the color screening constant to be determined 
by fitting the deuteron mass in this model. The quark 
delocalization in QDCSM is realized by allowing the sin- 
gle particle orbital wave function of QDCSM as a linear 
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combination of left and right Gaussian, the single parti- 
cle orbital wave functions in the ordinary quark cluster 
model, 
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The mixing parameter e(5) is not an adjusted one but 
determined variationally by the dynamics of the multi- 
quark system itself. This assumption allows the multi- 
quark system to choose its favorable configuration in the 
interacting process. It has been used to explain the cross- 
over transition between hadron phase and quark-gluon 
plasma phase jlfjj. The model parameters are fixed as 
follows: The u, d-quark mass difference is neglected and 
m u =nid is assumed to be exactly 1 /3 of the nucleon mass, 
namely m u —nid—m MeV. The it mass takes the exper- 
imental value. The A takes the same values as in Ref. Q , 
namely A=4.2 fm . The chiral coupling constant a c h 
is determined from the ttNN coupling constant as usual. 
The other parameters b, a Cl Vq, and a s are determined 
by fitting the nucleon and A masses and the stability 
of nucleon size b with the variation of quark mass m. 
All parameters used are listed in Table Q] In order to 
compare the intermediate-range attraction mechanism, 
the a meson exchange in ChQM and quark derealiza- 
tion and color screening in QDCSM, the same values of 
parameters: b, a s , a c h, m u , m^, A are used for these 
two models. Thus, these two models have exactly the 
same contributions from one-gluon-exchange and 7r ex- 
change. The only difference of the two models is coming 
from the short and intermediate-range part, a exchange 
for ChQM, quark derealization and color screening for 
QDCSM. To show the sensitivity of the QDCSM to the 
model parameters, the results of another set of model 
parameters (QDCSM2) is also reported. 



C. Quark delocalization model with hidden color 
channels coupling (QDCCM) 

This approach is focused on the hidden color chan- 
nel effect which has been ignored almost in all quark 
model calculations but certainly should exist in a de- 
scription based on the fundamental quark-gluon degree 
of freedom. In the lattice QCD calculation of NN inter- 
action 0] these hidden color channels should have been 
included implicitly. However their effect has not yet been 
separated. We assume a Hamiltonian which is the same 
as that of QDCSM except that the usual quadratic con- 



TABLE I: Parameters of three quark models discussed in this 
paper. 

ChQM QDCSM1 QDCSM2 QDCCM 



m u , d (MeV) 


313 


313 


313 


313 


&(fm) 


0.518 


0.518 


0.60 


0.518 


a c (MeVfm- 2 ) 


46.938 


56.755 


18.5 


56.755 


V (im 2 ) 


-1.297 


-0.5279 


-1.3598 


-0.5279 


Mfm" 2 ) 




0.45 


1.00 






0.485 


0.485 


0.996 


0.485 


m^MeV) 


138 


138 


138 


138 


Olch 


0.027 


0.027 


0.027 


0.027 


m CT (MeV) 


675 








Affm" 1 ) 


4.2 


4.2 


4.2 


4.2 



finement 



V C (r tJ ) 



-ka c Xi ■ Aj(r| + V ). 
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is used but with an additional multiplying factor k. For 
the color-singlet channels (two baryon clusters are in the 
color-singlet states), the factor k takes the value 1. For 
the hidden color channels, two recipes are used. Recipe 
1 (QDCCM1): For the coupling between hidden color 
channels and the color singlet channels, the factor k is 
taken as an adjustable parameter. All the other cases the 
factor k is kept 1. Recipe 2 (QDCCM2): The factor k is 
taken as adjustable parameter not only for color singlet- 
hidden color channels coupling but also for hidden color- 
hidden color channels. As for the single quark orbital 
wave function, the same form Eq.(3) as that of QDCSM 
is assumed. This model assumption is inspired by the 
lattice QCD calculation : The recent lattice QCD calcu- 
lations show that the interactions among quarks are gen- 
uinely multi-body interactions. The color dependent two 
body confinement interaction is consistent with the lat- 
tice QCD results only for two and three quark systems in 
color singlet states but inconsistent with the multi-body 
interaction obtained in lattice QCD for multi-quark sys- 
tems [l3| ■ So the direct extension of the color dependent 
two body confinement interaction from two- or three- 
quark system to multi-quark system as used in the most 
quark model calculations is questionable. The calculation 
based on the direct extension can not describe the NN 
scattering quantitatively well even after including hidden 
color channels coupling as shown in QDCCM0 might be 
an indication of this inadequacy. In fact, for multi-quark 
systems and color octet nucleons, quark pairs are not 
always in color antisymmetric state but also color sym- 
metric ones. The color factor A,- Xj will give rise to anti- 
confinement interaction for symmetric quark pairs [l8j . 
In QDCSM mentioned above, we used a color screening 
confinement interaction to model the effect of this multi- 
body confinement interaction obtained in lattice QCD. 
Here we study directly the effect of hidden color channel 
coupling to test if the phenomenological color screening 
confinement is an effective description of the hidden chan- 



4 



TABLE II: The channels used in NN scattering calculations 
and the factors k\, k 2 (for recipes 1,2) for each channel (1=1). 
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nel coupling. In order to simplify the numerical calcula- 
tion, a two body confinement interaction form Eq.(4) is 
still assumed but with an additional adjustable multiply- 
ing factor aimed to reflect the effect of the lattice QCD 
multi body confinement. At the same time, the model 
parameters are kept to the same as those of QDCSM1, 
except the color screening confinement form Eq.(2) is re- 
placed by the usual quadratic confinement form Eq.(4). 
This is aimed to let the effect of hidden color channel 
coupling stand out. 

III. THE RESULTS AND DISCUSSIONS 

We calculated the NN scattering phase shifts of dif- 
ferent partial waves (S, P, D, F, G, H and I waves) by 
three quark models mentioned above. To look for non- 
strange dibaryon resonances, a systematic calculation of 
NN scattering phase shifts with explicit coupling to iVA 
and AA channels is also done. The resonating-group 



TABLE III: The channels used in NN scattering calculations 
and the factors hi, k 2 (for recipes 1,2) for each channel (1=0). 
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method (RGM), described in more detail in 
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FIG 1. The phase shifts of NN 5-wave scattering. 

Ref.[l!|, is used to do the calculation. The experimental 
information used for the comparison is the partial-wave 
solution SP07 of NN scattering data. For QDCSM, 
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FIG. 2: The phase shifts of NN P wave scattering. 



the color screening parameter [i is fixed by deuteron prop- 
erties and no other parameters readjusted. For the third 
approach (QDCCM), the channels included in different 
partial waves are listed in Table II and III. The multiply- 
ing factors ki, k 2 are adjusted to fit the NN phase shifts 
of SP07. The calculated results for NN scattering phase 
shifts are presented in section A; deuteron properties are 
shown in section B and the discussions of the dibaryon 
resonances are given in section C. 



A. NN scattering phase shifts 

(1) S— waves: Fig. 1 shows the NN scattering phase 
shifts for 3 Si and 1 So partial waves. A perfect fit is ob- 
tained for both ChQM and QDCSM1 (QDCSM2 gives a 
little less attraction). The dominant contribution to the 
S-wave phase shift comes from the central part of the 
potentials. The agreement between two models means 
these two quark models give the same NN interaction, 
at least the same central part. For QDCCM, QDCCM1 
and QDCCM2 also give good descriptions of NN 3 Si and 
1 So scattering phase shifts by including the hidden color 
channels and adjusting the color confinement interaction 
strength, while with the usual color confinement interac- 
tion strength (k = 1), the model (QDCCMO) calculated 



phase shifts are far from the measured ones. 

(2) P— waves: Fig. 2 shows the NN scattering phase 
shifts of 1 P\, 3 P , 3 -Pi and 3 Pi partial waves. For 1 P 1 
and 3 Pi, ChQM and QDCCM gave an almost perfect de- 
scription of the experimental data. The 1 Pi phase shift 
is mainly determined by the central repulsion. The the- 
oretical phase shifts of QDCSM and QDCCMO are lower 
than experimental ones which show that these two mod- 
els give a too strong repulsion. For 3 P , QDCSM2 de- 
scribed the experimental data better than others. For 
QDCCM we do not have to adjust the color confine- 
ment interaction strength a c too much (fc = 1.1 for both 
QDCCM1 and QDCCM2), so both QDCCM1 and QD- 
CCM2, even QDCCMO, can fit the 3 P phase shifts rea- 
sonable well. ChQM and QDCSM1 give too strong at- 
traction. For 3 P 2 , QDCCM1 gives a perfect fit. ChQM, 
QDCSM and QDCCMO do not have enough attraction. 
Fig. 3 shows central, spin-orbit and tensor components of 
the 3 Pj phase shifts. Clearly, ChQM, QDCSM and QD- 
CCMO do not give strong enough attraction in the central 
and spin-orbit parts. In OBE, np, 7rw-exchange, which 
might have not been reproduced in the quark model cal- 
culations, are also needed to reproduce the P-wave phase 
shifts Q. 

(3) D— waves: Fig. 4 shows the NN scattering phase 
shifts of 3 D\, 3 D 2 , 3 D 3 and 1 D 2 partial waves. For 3 D\, 
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FIG. 3: Central, spin-orbit and tensor components of NN P 
wave scattering. 



all the models fit the experimental data well except QD- 
CCMO. For 3 D 2 , QDCSMl and QDCSM2 give a very 
good description of the experimental data (QDCSM2 is a 
little better) , ChQM gives too strong attraction. For QD- 
CCM, we find that we do not need adjust the color con- 
finement interaction strength for this channel, QDCCMO 
can fit the experimental scattering phase shifts. For 3 D 3 
and 1 D 2 , ChQM described the experimental data better 
than QDCSM. For QDCCM, both adjusting recipes can 
give a perfect fit to the experimental data. 

(4) F-wave: The calculated 1 F 3> 3 F 2 , 3 F 3 and 3 F 4 NN 
phase shifts are shown in Fig. 5. For F wave scattering, 
we find that QDCCMO already fit the experimental scat- 
tering phase shifts reasonably so we did not fine tune the 
confinement strength . All the models give a good de- 
scription of the experimental data reasonably in the low 
energy region (E c . m . < 100 MeV). Above 100 MeV, the 
model predictions deviate more or less from the experi- 
mental data. For 3 F 2 QDCSM2 gave much better fit to 
the experimental data than QDCSMl, ChQM and QD- 
CCMO. For l F 3 , QDCSMl, QDCSM2 and QDCCMO all 
give better fit to the experimental data than ChQM, es- 
pecially at higher energy. However, for 3 i<3, ChQM is 
closer to the experimental data than other models. For 
3 F 4 , a perfect fit is obtained for QDCCMO, ChQM has a 
little too strong attraction at high energy and QDCSM 
gives a too weak attraction. 

(5) G-wave: The NN phase shifts of 3 G 3 , 1 G A , 3 G 4 and 
3 G5 are shown in Fig. 6. All the models can describe the 
experimental data. We do not have to adjust the color 
confinement interaction strength for QDCCM here. 

(6) //-wave: Fig. 7 shows the calculated 1 H 5 , 3 H 4 , 3 H 5 
and 3 Hq NN phase shifts. For //-wave phase shifts, all 
models fit to the experimental data equally well. We also 
find that we do not have to adjust the color confinement 
interaction strength for QDCCM here. 



(7) /-wave: The calculated 3 /5, 



3 I 6 and 3 I 7 NN 



phase shifts are shown in Fig. 8. For /-wave phase shifts, 



all the models give almost the same results and fit the ex- 
perimental data well. Again the color confinement inter- 
action strength for QDCCM do not need to be adjusted 
here. 

For high L partial waves, the long range it exchange 
dominates the interaction. Three quark models have the 
same tt exchange and therefore they give almost the same 
results for L > 3 and we do not have to adjust the mul- 
tiplying factor for the QDCCM for these high L partial 
wave. 

These numerical results (Figs. 1-8) show that by includ- 
ing the hidden color channels and adjusting the color 
confinement interaction strength, both adjusting recipes 
can fit the JVJV scattering phase shifts well. From the 
calculated S, P, ZD-wave phase shifts of AW scattering in 
QDCCMO, we can see that the attraction is always in- 
adequate because of the appearance of anti-confinement 
interaction of symmetric quark pairs. By increasing the 
strength of confinement, the attraction coming from the 
confinement interaction is strengthened, QDCCM1 and 
QDCCM2 can give a good description of the experimen- 
tal data. We take these results as an indication that the 
short and intermediate range AW interaction is caused 
by the nucleon internal structure and its distortion both 
in orbital and color spaces in the interacting process. 
These are quite the same as the atomic internal struc- 
ture and its distortion in orbital space which give rise 
to the molecular covalent bond. The Anderson's con- 
jecture Q is verified here. The phenomenological color 
screening confinement might be an effective description 
of the hidden color channel coupling. The phenomenolog- 
ical a meson exchange used in OBE and ChQM might be 
an effective description of the more complicated nucleon 
distortion in the NN interaction process as described in 
QDCSM and QDCCM. This mechanism also gives a nat- 
ural explanation why does the AW interaction between 
two color singlet nucleons is so similar to the molecular 
interaction between two charge neutral atoms except the 
energy and length scale difference. 



B. Deuteron 

All these three models are used to calculate the proper- 
ties of deuteron, the results are shown in Table IV. Both 
ChQM and QDCSM give a good description of deuteron. 
For QDCSM, by adjusting the color screening parame- 
ter, the same results for deuteron can be obtained for 
different baryon size b. Because of the large separation 
between the proton and neutron in the deuteron, the 
properties of deuteron mainly reflect the long-range part 
of the nuclear force. The same 7r-exchange used in the 
two models assure the properties of deuteron be fitted 
equally well. However ir exchange alone can not provide 
strong enough intermediate-range attraction to make the 
deuteron bound. In ChQM, it is the phenomenologi- 
cal a meson exchange which provides the intermediate 
range attraction. In QDCSM it is the quark delocal- 
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FIG. 4: The phase shifts of NN D wave scattering. 
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FIG. 5: The phase shifts of NN F wave scattering. 
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FIG. 6: The phase shifts of NN G wave scattering. 
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FIG. 7: The phase shifts of AW H wave scattering. 
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FIG. 8: The phase shifts of NN I wave scattering. 



ization and color screening which provide the intermedi- 
ate range attraction. The fact that both models fit the 
deuteron properties well verifies once more that the two 
intermediate range attraction mechanisms used in these 
two models are equivalent. 

Table IV shows that the binding energy and the 
D-wave component of deuteron can be reproduced (we 
didn't fine tune the strength of color confinement to get 
a better fitting). However the root mean square radius is 
too small in comparison to experimental value. This may 
indicate that QDCCM with the parameters giving in Ta- 
ble I gives rise to an NN scattering phase shift equivalent 
potential but a little too strong attraction in the short 
range region, which tightens up the deuteron. 

TABLE IV: The properties of deuteron. 

ChQM QDCSM1 QDCSM2 QDCCM1 QDCCM2 

B (MeV) 
vr^(fm) 
Pd(%) 



C. Dibaryon resonances in NN scattering 

In this part, we show the results of a systematic search 
for the possible non-strange dibaryon candidates by three 
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quark models mentioned above. 

The previous calculations [2l[ show that there are four 
possible dibaryons in the quark model calculations, NA 
state with I J = 12, AA states with IJ = 01, 10, 03. Here 
the QDCCM is applied to recalculate these states. All of 
these dibaryon states are allowed to decay via the NN 
channels. In other words, these dibaryon states appear 
as resonance states in the 7V7V scattering process. So we 
calculate the NN scattering phase shifts by including all 
the possible channel couplings. The results are shown in 
Fig. 9. 

(1) J = 0, J = 1: The 3 Si energies of single A— A chan- 
nel calculation are lower than the corresponding thresh- 
old 100-350 MeV in ChQM, QDCCM and QDCSM. The 
coupling to the 3 S^ N channel has an unexpectedly large 
effect, pushing up the energy of 3 S^~ A state ~ 300 MeV, 
so that only in QDCSM2 it becomes a resonance at 2408 
MeV. This very large mass shift is caused by the central 
interaction and the presence of a lower-mass state, the 
deuteron, in the admixed 3 Si N channel. Mixing with 
other channels listed in Table III, the resonance mass is 
pushed down a little bit, to 2393 MeV. In ChQM and 
QDCCM, the 3 Si energies in the single A — A channel 
calculation are 100 MeV or more higher than that in QD- 
CSM2. The additional large mass shift caused by the cou- 
pling to the NN channel then pushes the state above the 
AA threshold. So no resonance appears in other models 
except QDCSM2. The phase shifts of 3 S NN are shown 
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FIG. 9: The phase shifts of NN S wave and D wave scattering to energies beyond the AA or JVA threshold. 



in the up left corner of Fig. 9, where the phase shifts for 
100 < E c . m . < 400 MeV from ChQM, QDCSM1, QD- 
CSM2, QDCCM1 and QDCCM2 agree with each other, 
as already pointed out in section A. The phase shifts 
of 3 S^ N rises through it/ '2 at a resonance mass only in 
QDCSM2. So, is there an IJ — 01 AA resonance state 
with resonance mass 2393 MeV in the NN 3 S± scattering 
channel is not sure. 

(2) I = 0, J = 3: The single A — A channel calculation 
shows that the state 7 5 AA is a bound state in all models 
used here. The coupling to the 3 D^ N channels causes 
this bound state change into an elastic resonance. The 
resonance mass shift, which is caused by the tensor inter- 
action, is not large ~ 3 MeV. The calculation shows that 
the mass shift is always dominated by the NN scattering 
states below the bound-state rather than those above it. 
Coupling to other channels listed in Table III which are 
above the 7 S AA bound state, the resonance is pushed 
down as expected. The calculated 3 D^ N phase shifts, 
shown in the up right corner of Fig. 9, rise through 
it/2 at the resonance masses in all models. But quan- 
titatively the resonance masses are different in different 
models. The resonance mass in QDCSM1 is about 60 
MeV lower than that in ChQM, and the QDCSM2 al- 
ways has the lowest mass. For QDCCM, the resonance 
mass is 2443 MeV in QDCCM0, 2298 MeV in QDCCM1 
and 2156 MeV in QDCCM2. This resonance (I J = 03 
AA) is a promising candidate for the observed isoscalar 



ABC structure seen more clearly in the pn — > dmr pro- 
duction cross section at 2.36 Gev in the recent report by 
the CELSIUS- WAS A Collaboration 

(3) I = 1,J = 0: The X S AA state is qualitatively 
similar to the 3 S' AA state, since they are just different 
spin-isospin states of the same quark system with the 
same relative orbital angular momenta. The calculated 
phase shifts, shown in the down left corner of Fig. 9, 
show that the resonance survives only in QDCSM2 after 
the channel coupling. The situation is almost the same 
as the 3 S^ N state. 

(4) I = 1, J = 2: The phase shifts of NN scattering 
are shown in the down right corner of Fig. 9. From the 
curves, we find that a resonance appears in QDCSM2, 
QDCCM1 and QDCCM2. The resonance masses are: 
2168 MeV in QDCSM2, 2144 MeV in QDCCM1 and 2130 
MeV in QDCCM2. For ChQM and QDCSM1, only a 
prominent cusp appears at the NA threshold. Neverthe- 
less, the state might correspond to the resonance looping 
in the Argand diagram of the 1 D 2 pp-partial wave [23j. 

For odd-parity NN states, resonance poles are found 
for the isovector odd-parity NN partial waves 3 P 2 , 3 F 2 
and 3 F$ 24). These empirical resonance-like solutions 
reproduce the empirical Argand loopings of the partial 
wave solutions, but many studies in the past [25j have 
not resolved the question of whether these Argand loop- 
ings represent real dibaryon resonances. In our quark 
model calculation, we have not found any resonance at- 
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tributable to an AA or AA bound state in the odd-parity 
NN states. 



IV. SUMMARY 

By including the hidden color channels and varying 
the strength of the color confinement potential between 
color-singlet channels and hidden color channels and/or 
hidden color channels and hidden color channels, a phe- 
nomenological quark model for baryon-baryon interac- 
tion is constructed. The model achieves a good descrip- 
tion of S-, P-,D-,F-,G-,H-, /-partial wave phase shifts of 
NN scattering as good as other quark models. It also 
reproduces the binding energy and D-wave component of 
deuteron but a little too small root mean square radius. 
Applying the model to dibaryon search, similar results 
with QDCSM and ChQM are obtained. The results show 
that the hidden color channels are important for the NN 
intermediate range attraction. The lattice QCD calcu- 
lations obtained the string like multi-body confinement 
interaction in the multi-quark system [l7| . It is equiva- 
lent to the two body confinement Eq.(4) with k = 1 for 
a color singlet nucleon with three quarks. Oka extended 
the string-flip model to six-quark system and obtained 
a reasonable description of NN interaction [llj which 
might be viewed as a modeling of the lattice QCD string 



like multi-body confinement. QDCCM fits the NN scat- 
tering data better and we suspect it might be another 
modeling of the lattice QCD multi-body confinement. 

Certainly one would expect to directly use the string 
like multi-body interaction obtained in lattice QCD to 
calculate the NN interaction. However, it is not only be- 
cause of the huge numerical task but also because there is 
no any information about the transition interaction be- 
tween different string structure which hindered this ap- 
proach. 

Nuclear force is an old topic, it has been studied 
over 70 years and a large amount of experimental data 
has been accumulated. Although there are several ap- 
proaches which can give almost perfect description of the 
experimental data, the mechanism for the intermediate- 
range attraction is still an open question. Lattice QCD 
achieved a qualitative description of the A^A" interaction 
already and it will finally achieves a quantitative descrip- 
tion. But based on present lattice QCD technique it 
can not reveal the physical mechanism, for example to 
distinguish the phenomenological a meson exchange and 
the nucleon distortion similar to molecular covalent bond 
mechanism for the intermediate range attraction. One 
has to develop the non-perturbative continuous QCD 
field theory method as well as non-perturbative QCD 
model to explore the NN interaction. 
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